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Method and arrangement for 
medical X-ray imaging 

BACKGROUND OF THE INVENTION 

Three-dimensional X-ray imaging is based on taking 
several 1-D or 2-D projection images of a 3-D body 
from different directions. If 1-D projection images are 
available from all around a 2-D slice of the body with 
dense angular sampling,, the inner structure of the slice 
can be determined. This is known as Computerized 
Tomography (CT) imaging technology, which is widely 
used in medicine today. A crucial part of CT technology 
is the reconstruction algorithm taking the X-ray images 
as argument and producing a voxel representation of 
the 3-D body. 

In many practical cases X-ray projection images are 
available only from a limited angle of view. A collection 
of X-ray images of a 3-D body is called sparse 
projection data if (a) the images are taken from a 
limited angle of view or (b) there are only a small 
number of images. Sparse projection data does not 
contain sufficient information to completely describe the 
3-D body. 

However, some a priori information about the body is 
typically available without X-ray imaging. Combining 
this information with sparse projection data enables 
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more reliable 3-D reconstruction than is possible by 
using only the projection data. 

Traditional reconstruction algorithms such as filtered 
backprojection (FBP), Fourier reconstruction (FR) or 
algebraic reconstruction technique (ART) do not give 
satisfactory reconstructions from sparse projection 
data. Reasons for this include requirement for dense 
full-angle sampling of data, difficulty to use a priori 
information, for example nonnegativity of the X-ray 
attenuation coefficient, and poor robustness against 
measurement noise. For example the FBP method relies 
on summing up noise elements with fine sampling, 
leading to unnecessarily high radiation dose. 

The traditional setting of the image reconstruction 
problem in tomography assumes that the imaging 
geometry is known accurately. The prior art requires 
either reference feature points, with which the unknown 
part of the imaging geometry (e.g. motion of the 
imaging device or the object) is computed, or precise 
prior knowledge of the imaging geometry. Anyway the 
prior art does not know effective methods to determine 
the imaging geometry like the motion of the imaging 
device and/or the object, though all parts of the 
imaging geometry are rarely precisely known in 
advance. 
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BRIEF DESCRIPTION OF THE INVENTION 

The aim of the invention is to overcome the problems 
met in 3-D reconstruction of the object that occur when 
using traditional reconstruction algorithms with sparse 
projection data. This is achieved by a method for 
producing three-dimensional information of an object in, 
medical X-ray imaging. The object is X-radiated from at 
least two different directions and the said X-radiation is 
detected to form projection data of the object. The 
object is modelled mathematically utilizing the 
projection data to solve the imaging geometry and/or 
the motion of the object. Said projection data and said 
mathematical modelling of the object are utilized in 
Bayesian inversion based on Bayes' formula 



to produce three-dimensional information of the object, ' 
the prior distribution Ppr(e) representing the prior 
knowledge of the imaging geometry and/or the motion 
of the object, the prior distribution p pr (x) representing 
mathematical modelling of the object, x representing 
the object image vector, which comprises values of the 
X-ray attenuation coefficient inside the object, 9 
representing the parameter vector of the imaging 
geometry and/or the motion of the object, m 
representing projection data, the likelihood distribution 
p(m|x,6 ) representing the X-radiation attenuation 
model between the object image vector x, geometry 
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parameter vector B and projection data m, p(m) being 
a normalization constant and the posteriori distribution 
p(x,6|m) representing the three-dimensional 
information of the object and the imaging geometry 
including the motion of the object. 

The invention also relates to a medical X-ray device 
arrangement for producing three-dimensional 
information of an object (4) in a medical X-ray imaging, 
The medical X-ray device arrangement comprises: 

- an X-ray source for X-radiating the object from at 
least two different directions 

- a detector for detecting the X-radiation to form 
projection data of the object 

means for modelling the object mathematically 
utilizing the projection data to solve the imaging 
geometry and/or the motion of the object, where the 
solving concerns either some or all parts of the imaging 
geometry and/or the motion of the object. 

and the medical X-ray device arrangement compnses: 

- means for utilizing said projection data and said 
mathematical modelling of the object in Bayesian 
inversion based on Bayes' formula 



p(x 7 &\m) = 

pirn) 



to produce three-dimensional information of the object, 
the prior distribution p pr (0) representing the prior 
knowledge of the the imaging geometry and/or the 
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motion of the object (4), the prior distribution p pr (x) 
representing mathematical modelling of the object, x 
representing the object image vector, which comprises 
values of the X-ray attenuation coefficient inside the 
object, e representing the parameter vector of the 
imaging geometry and/or the motion of the object (4), m 
representing projection data, the likelihood distribution 
p(m[x,e ) representing the X-radiation attenuation 
model between the object image vector x. geometry 
parameter vector 9 and projection data m, p(m) being 
a normalization constant and the posteriori distribution 
p(x,e|m) representing the three-dimensional 
information of the object (4) and the imaging geometry 
Including the motion of the object- 
In this invention, the imaging geometry and/or the 
motion of the object are solved with the image 
reconstruction simultaneously or having different 
timings. This is based on that projection data itself 
simultaneously contains sufficient amount of information 
about the imaging geometry and the 3D reconstruction. 
Neither reference feature points nor prior knowledge of 
the motion are needed. 

The invention is also based on that biological issues 
have that kind of statistical a priori information that this 
information can be successfully utilized with Bayesian 
inversion in medical x-ray imaging. The suitable a priori 
information makes possible to model the biological 
tissue mathematically accurately enough and 
independently of X-ray imaging. From biological tissue 
it ie possible to compile qualitative structural information 
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which makes it possible to utilize the Bayesian method 
successfully to solve the problems in medical three- 
dimensional x-ray imaging. There is certain regularity in 
biological tissues and this regularity is useful especially 
with the Bayesian method. 

For example 1 0 x-ray images are taken from breasts of 
different persons. From these x-ray images it is noticed 
that there is much similarity in the statistical structure of 
the breasts between different people. In other words 
biological tissues and x-ray images taken from the 
biological tissues have similar or almost similar 
statistical structure between different persons. 

With Bayesian inversion it is possible to utilize a priori 
information efficiently in 3-D reconstruction from sparse 
projection data, because the suitable a priori 
information from the biological tissue makes possible to 
model the biological tissue mathematically accurately 
enough and independently of X-ray imaging. In 
addition, a priori information of the imaging geometry, 
such as approximate positions of the radiation source 
and the detector, can be similarly taken into account in 
a sound, statistical way. Any collection of projection 
data can be used in the reconstruction. Application- 
dependent a priori knowledge can be used to regularize 
the ill-posed reconstruction problem. 

This invention improves the quality of 3-D 
reconstructions over traditional methods. In addition, 
the number of radiographs can be minimized without 
compromising quality of the reconstruction. This is very 
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important in medical applications since the X-ray dose 
of the patient can be lowered. 



BRIEF DESCRIPTION OF THE FIGURES 

Figurel shows a simple example of X-ray imaging. 

Figures 2A -2B show a difference between global 
tomography and local tomography. 

Figures 3A - 3B show different types of sparse 

projection data. Every black dot represents 
a location of the X-ray source for taking 
one projection image. 

Figures 4A - 4B show examples of parts of boundaries 
and cracks of an object visible and invisible 
without a priori information. 

Figure 5 illustrates "pencil beam" X-ray attenuation 
model. 

Figure 6 is basic flow chart of the method according 
to the invention. 

Figure 7 is an intraoral X-ray device arrangement 
presenting one preferred embodiment of 
the invention. 
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Figure 8 shows a measurement geometry of dental 
limited-angle tomography with a digital 
intraoral sensor. 



Figure 9 shows results of 3D reconstruction of head 
phantom in the first preferred embodiment. 



DETAILED DESCRIPTION OF THE INVENTION 

In conventional tomography, the interior of an object is 
reconstructed from tomographic projections such as X- 
ray images. All the current image reconstruction 
methods assume that projection geometry of the 
imaging device is either known or solved in advance by 
using fiducial or non-fiducial feature points in the 
images. With this invention a novel approach is 
presented where the imaging geometry is solved with 

the reconstruction problem while no correspondence 
information is needed. This approach is a direct 
application of Bayesian inversion theory and produces 
the maximum likelihood or maximum a posteriori 
estimates for the imaging geometry like motion 
parameters under the selected noise and prior 
distributions. This invention can be utilized 
assuming either ID or ZD projections. To utilize the 
invention with ID projections is also often a sufficient 
assumption since the reconstruction problem is 
frequently performed as a stack of 2D reconstruction 
problems for computational reasons. However, the 
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method can be also directly extended to the general 3D 
case. 

In practical imaging situations X-ray images are not 
always available from all around the body. The body 
might be visible only from certain directions due to 
imaging geometry. For example this is the case in 3-D 
mammography with the breast compressed against the 
detector, or in intraoral dental imaging with the 
detector inside the patient's mouth. This situation is 
called limited-angle tomography. Also, the region of 
interest might be surrounded by tissue that need not be 
imaged, like in extrapral dental imaging. This situation 
is called local tomography. In addition the number of 
radiographs should be minimized in medical 
applications for reducing the X-ray dose of the patient. 

In the preferred embodiments of the invention Bayesian 
inversion algorithms are used to create a new type of 
3-D medical X-ray imaging. It is intermediate between a 
projection radiograph and a full computed tomography 
scan. Two steps are needed to perform successfully 
such imaging: In step one, the doctor (a) collects all 
the prior information he or she has on the tissue of 
interest, and (b) takes the minimum number of 
radiographs containing the desired 3D information. In 
step two, a tomographic algorithm based on Bayesian 
inversion is used for 3D reconstruction. The algorithm 
contains two separate mathematical models. First, all a 
priori information (i.e, information that is independent 
of X-ray imaging) is used to model the unknown object 



mathematically. The model is put into the form of prior 
distribution in Bayes' formula- In addition, prior 
information of the imaging geometry parameters can be- 
modelled. Second, the measurement is modelled 
mathematically. This involves the geometry of the 
imaging device, position of the detector and x-ray 
source during the exposure of each projection image 
and a model for X-ray attenuation inside tissue. The 
mathematical model of the measurement is put into the 
form of likelihood distribution in Bayes' formula. 

In figure l is shown a simple example of x-ray imaging, 
where an X-ray source 2 is placed on one side of an 
object 4 under imaging. Radiation passes through the 
object and is detected by a detector 6 on the other 
side. The X-ray source is for example a X-ray source of 
an Intraoral X-ray source of a dentist, of a dental 
panoramic X-ray device, of a surgical C-arm X-ray 
device, of a mammography device or of any other 
medical X-ray device and the detector 6 is a detector of 
some of those devices, usually the detector 6 is a 
digital sensor that can be thought of as a 2-D array 
almost pointlike detectors. 
The 3-D body under imaging is modelled by 
nonnegative X-ray attenuation coefficient. The value 
gives the relative intensity loss of the X-ray travelling 
within a small distance dr; 

dI < w 
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(1) 



The X-radiation has initial intensity / 0 when entering the 
object 4 and a smaller intensity I x when exiting the 
object. The following equation shows the attenuation 
law: 

f f I'ir) 

I x(r)dr = — I 77Y<ir = logi" 0 — l°g-7i 
Jl Jl I K r ) 

(2) 

where initial intensity / 0 is known by calibration, and 

intensity after object /j is known from the 
corresponding point value in a projection image. Thus 
the measured data is the integral of x along the line L 

In the above model it is not taken into account (a) 
scattering phenomena resulting in X-rays changing their 
direction, or (b) the dependency of attenuation on the 
X-ray spectrum resulting in low-energy photons being 
more easily attenuated than high-energy ones. Effect 
(b) causes errors in measurements and is sometimes 
referred to as beam hardening. More detailed models 
can be utilized in Bayesian inversion to overcome (a) 
and (b). 

In medical imaging the geometrical arrangements of 
the X-ray source and digital sensor vary according to 
the diagnostic task and equipment. Figures 2A-2B and 3 
illustrate the types of tomographic data resulting from 
different imaging situations. For clarity, here are 



presented two-dimensional examples; similar situations 
can be considered in 3-D. 

In figure 2A-2B are shown two cases according to 
whether the whole object 4 is fully visible In each 
projection or not. These cases are called global 
tomography and local tomography, respectively. In 
figure 2B is radiated only ROI (Region Of Interest) 8- 

The types of data described in figures 3A - 3B cover a 
large range of specific imaging tasks. The choice of 
data collection dictates what kind of features and 
details it is easiest to reconstruct reliably from the data. • 
Thus the choice of projection views must be made 
according to the diagnostic task at hand. In traditional 
CT imaging, projections are taken from all around the 
object 4. In the preferred embodiments of the invention 
radiation dose is lowered by sampling the angular 
variable more sparsely. Figure 3A presents more dense 
angular data with full angle and limited angle 
measurements. Figure 3B presents more sparse angular 
data, also with full and limited angles of measurement. 

In figures 4A and 4B there are presented examples of 
parts of boundary and cracks that are visible or 
indetectable in reconstruction without a priori 
information. Figure 4A there is an object 4 under 
imaging with edge on the surface of the object- The 
edge is detectable on the leftmost position where the 
edge is more parallel to the direction of the X-rays. In 
the figure 4B there is a crack inside of the object 4. The 
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crack is detectable in the leftmost position parallel to 
the angle of measurement.' 

The main idea in Bayesian inversion method is to 
consider the inverse problem as a problem of statistical 
inference. All variables are redefined to be random 
variables. The randomness reflects uncertainty of their 
actual values and the degree of uncertainty is coded in 
the probability distributions of these random variables. 

Typically, the mode) for X-ray tomography problem is 
represented by the formula 

(3) m~Ax+e, 

where m is a vector of projection measurements (all 
the projetion images represented as one vector), x is a 
vector that represents the object densities (voxel 
values) at given positions in 3D space, A is a 
transformation matrix that characterises how 
the projections are constructed from the 3D object, 
hence contains the imaging geometry, and e is a 
random, noise vector. The variables m, x and e are 
random variables and the matrix A comes from the 
pencil beam model for the X-ray attenuation: A contains 
the weights describing how much each voxel 
contributes to the X-ray attenuation along the X-ray 
paths. In this way the integral in formula (2) is 
approximated. This is illustrated by an example in 
figure 5, 
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In the traditional setting, the objective is to reconstruct 
the structure (density) of the object given the 
projection images and the imaging geometry- As an 
extension to the traditional setting, in this invetion both 
the reconstruction and the imaging geometry are 
solved simultaneously. To emphasise that, we write 
the model (3) in the form 

m = Aq x + e, 

where we just explicitly say that A depends on the 
geometry parameters e that define the unknown part of 
the imaging geometry and motion of the object. Hence, 
the unknowns in the equation are both 0 and x. 

Simultaneous solving of both 8 and x is possible, 
because, the model (3) constrains the imaging 
geometry sufficiently- This is due that, even though 
there are more unknowns than equations, the 
equations are not independent and, if e.g. the motion 
parameters are incorrect, the equations are also 
inconsistent. Moreover, if the geometry parameters 8 
were incorrectly fixed, in general, there is no 
reconstruction x that wouid produce the observed 
images m. Therefore, excluding certain special cases, 
the geometry parameters are completely determined 
from the noise model, though the image object x is not : 
uniquely (additional prior information is required), as 
far as the imaged object is sufficiently complex. 
Therefore any feasible reconstruction x (in the 
likelihood sense) can be used to compute 8, after which 
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the reconstruction x can be updated to accord best with 
the prior information. 

Assume now that the object image vector x, geometry 
parameter vector 0, and the noise are jointly 
independent random variables. The conditional 
probability distribution of x given the measurement m is. 
given by Bayes' formula in the form 



pi_x,e\ni) = 



p pr (e)p pr (x)p(m\x,0) 



p{tri) 
(4) 

where p(m) is a normalization constant. The density 
P(x,e|m) is called the posteriori distribution of x and 0. 
p(m|x,9) is a likelihood distribution representing the X- 
radiation attenuation model between the object image 
vector x, geometry parameter vector 9, and projection 
data m. 



The density P pr ( x ) , called the prior distribution of x, 

is designed to contain all possible information available 
of the actual object 4 independently of X-ray imaging. 

The density P P X&) is the prior distribution of 9, and it 

models the available information of the imaging 
geometry including the motion of the object. If we do 
not have any such prior information of the parameters 
9, we may use the non-informative prior 

P P r(&) = constant . it is crucial, in contrast to several 
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classical regularization methods, that the choice of the 
prior distributions should not be based on the 
projection data, me rule of thumb in the design of prior 
distributions is that typical image vectors (say, of some 
existing library) should have high prior probability while 
atypical or impossible ones should have low or 
negligible probability. 

In the framework of the Bayesian inversion theory, the 
posterior distribution in formula (4) represents the 
complete solution of the 3D reconstruction problem as 
well as of the imaging geometry estimation problem. 

To produce an image of the object 4, and similarly an 
estimate for the geometry parameters, based on the 
posterior distribution, several alternatives exist. The 
most common ones are the maximum a posteriori 
estimator (MAP) and conditional mean estimator (CM). 
Let y = (x,6) i.e. y fs the vector composed of both the 
object image and geometry parameter vectors. The 
estimators are defined by the formulas 

P(3W I "0 = max pdy \ rri) 
(5) 

where the maximum on the right hand side is taken 
over all y, and 
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(6) 

Finding the MAP estimator is an optimization problem 
while finding the CM estimator is a problem of 
integration. 

In figure 6 is presented a basic flow chart of the 
method according to the invention. In method step 600 
the object is modelled mathematically Independently of 
X-ray imaging. In method step 602 the object is X- 
radiated from at least two different directions. In 
method step 604 the said X-radiation is detected to 
form projection data of the object. In method step 606 
said projection data and said mathematical modelling of 
the object are utilized in Bayesian inversion based on 
Bayes formula to produce three-dimensional 
information of the object. The method step 600 is also 
possible to perform after method step 602 or 604, 

In the first preferred embodiment of the invention is 
presented an application to dental radiology. 

X-ray projection images are conventionally used in 
dental radiology. However, certain diagnostic tasks 
require more precise knowledge of the 3D structure of 
tissue than is available in two-dimensional radiographs. 
Such tasks include implant planning and detection of 
bone loss between tooth roots. 
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In figure 7 is presented an intraoral X-ray device 5 
arrangement presenting the first preferred embodiment 
of the invention. It is important to notice that this is 
only an example of the medical X-ray device 5 
arrangement where the invention is possible to be 
utilized. 

The medical x-ray device 5 in the preferred 
embodiments of the invention is for example a dental 
panoramic X-ray device, a surgical C-arm X-ray device 
or a mammography device. 

In the first preferred embodiment of the invention the 
articulated arm arrangement 3 moves the X-ray source 
2 to the right position. The X-radiation begins by 
pressing the exposure button 12. The X-ray source 2 X- 
radiates the object 4, which is for example teeth of a 
patient- The detector 6 detects the X-radiation, The 
image information which is got by detecting the X- 
radiation is sent by communication link 16 to the 
computer 14. The computer comprises the software 
means 15 to process the image information according 
to the invention. There can be more than one computer 
14 and also the software means 15 can situate in more 
than one computer 14. For example the first computer 
14 is computer which is used in x-ray imaging. The 
second computer 14 is computer which is used in 
processing the image information according to the 
invention. It is possible to have the second computer 
14 far away from the actual medical x-ray device 5. For 
simplicity in figure 7 is shown only one computer 14. 
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In the first preferred embodiment (Fig. 7) of the 
invention, the dentist's X-ray equipment is used for 
taking a set of 2D projection images that are used as 
input for Bayesian 3D reconstruction algorithm. Such 
equipment includes an intraoral X-ray unit and a digital 
intraoral sensor. 

Benefits of this approach over conventional CT scan are 

- low cost and convenient usage, 

- high resolution of the projection images, 

- possibility to take as few radiographs as is needed 
to capture the relevant 3D information, minimizing 
radiation dose to the patient, 

- the possibility to choose imaging directions so that 
the X-rays do not pass through the whole head but 
only the interesting tissue, further reducing dose. 

In the preferred embodiments of the invention, 
Bayesian inversion is used for the 3D reconstruction- 
Input data for the algorithm is the set of projection 
images and the following a priori knowledge: 

(a) Dental tissue consists of few approximately 
homogeneous regions with well-defined boundaries. 

(b) Dental tissue can only attenuate X-radiation, not 
intensify it. 

The example of the detector 6 used In the first 
preferred embodiment of the invention is based on 
charge coupled device (CCD) technology and has 
dynamic range of 4096 gray levels. The size of the 
active imaging area is 34 mm * 26 mm and the 
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resolution is 872 * 664 pixels- After exposure, each 
pixel contains an integer which is proportional to the 
number of X-ray quanta that hit the pixel's area. 

Alternative detectors Include any other digital intraoral 
sensor, digitized X-ray film, or any intraoral sensing 
device converting detected X-ray photons to a digital 
image. 

In the measurement geometry the focal point of the X- 
ray source 2 moves horizontally on a circle in the plane 
10, with center at the detector 6, see Figure 8. in 
Figure 8, the detector 6 corresponds to one row of CCD 
pixels in the CCD sensor. 

In figure 8 is presented an example of taking seven 
radiographs of some teeth of a patient. The total 
opening angle of the projection views is 55 degrees in 
this example. Figure 8 represents the xy plane, and the 
z-coordinate axis is thought to be perpendicular to the 
plane of the paper. 

In the preferred embodiments, mathematical modelling 
of the object 4 can be done with a prior distribution of 

the form in formula (7): 

PprO) = exp(-oj ^2 UwW) 

N 

(7) 

where the sum is taken over a collection of 3D 
neighbourhoods N and the value U N (x) depends only 
on the values of voxels belonging to the neighborhood 
N, and a is a positive regularization parameter used to . 
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tune the width of the prior distribution. The 
neighborhood N is typically a set of voxels whose center 
points are closer to each other than some predefined 
maximum distance. 

It is to be noted that formula (7) does not define a 
probability distribution since it does not Integrate to 1. 
Instead, the integral diverges. However, when used in 
the Bayes formula (4) with a suitable likelihood 
distribution, the posterior distribution is integrable. The 
same remark concerns the priors derived from (7), that 
is, formulae (9), (10) and (11). 

In the preferred embodiments of the invention the 
three-dimensional problem can be reduced to a stack of 
two-dimensional problems each corresponding to a 
plane determined by a constant value of z. Here figure 
_8 represents exactly the situation in the xy plane, i.e. 
the plane z=0. Each row in the detector corresponds to 
one such 2D problem. This approach leads to small 
approximation error because the X-ray source does not 
move in the correct plane for nonzero z, and this error 
is neglected. 

Next the modeling of a 2D tomographic problem is 
explained. Let 

(8 ) ma) = Ax(j) + e(j) ' 

denote the jth 2D tomographic problem- Here the 
vector m(j) contains the readings on jth row from each 
of the seven radiographs. The vector x(j) is the jth slice 
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of the 3-D representation x of the object 4 under 
imaging. The matrix A is the pencil beam model matrix 
for the 2-D tomography problem. Now x(j) is a 2D array 
of pixels. The pixels are denoted by x(j)[k,q], where 
k=l,2,...,K is row index and q=l,2,...,Q is column index. 

The mathematical modelling of the object 4, i.e. 
Incorporation of prior information, is next explained for 
the 2D slice. Define 

P^ a) ) =exp(-a^X3C7(iNT(o;^[fc > c7]))) 

k=l 5=1 

(9) 

where the value U(N(x(j)[k,q])) depends only on the 
values of pixels belonging to the neighborhood N of the 
pixel x(j)[k,q] and a is a positive regularization 
parameter used to tune the width of the prior 
distribution. TTie choice of the neighborhoods 
N(x(j)[k,q]) is arbitrary, but a typical choice is the set of 
pixels whose centerpoints are closer to the centerpoint 
of x(j)[k,q] than some predefined distance. 

In the preferred embodiment, the neighborhoods are 
chosen to consist of two adjacent pixels. Further, the 
function U is chosen to be a power of 
absolute value of the difference between adjacent 
pixels. These choices lead to formula (10): 

Ppr (o^) = exp ( - a( J2 \* U) &< d - + L «r] |* + 

+££i« w iM]-« ( ' ) [M+*ir 
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(10) 

where 5 is a positive real number and a is a 
regularization parameter used to tune the width of the . 
prior distribution. 

In the equation (11) is presented non-normalized total 
variation (TV) distribution, when s=l in equation (10). 
This is the case according to the preferred embodiment 
of the invention, because when s=l the prior model 
becomes such that it gives high probability for objects 
consisting of a few areas of different attenuation with 
well-defined boundaries. 



prTV(*(f>) - «p ( - 52 £ !•«[*, q) - + 1,?] | 4- 



(11) 

X-radiation can only attenuate, not strengthen, inside 
tissue- This leads to positivlty prior pos defined by 

(12) pos(x(j)) ^ 1 if all pixels of x(j) are positive, 0 
otherwise. 

The explanation of modeling of one 2D tomographic 
problem is now complete. 

The stack of 2D problems are connected to each other 
by demanding that consecutive 2D slices x(j) and x(j-l) 
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should not be very different. Mathematically this is 
expressed as follows: 

(13) pr3D(x(j)) = exp(- y ZZ I x(j)[k,q] - x(j-l)[k,q] 
I), 

where the sums are taken over k=l,...,K and q=l,...,Q 
and y > P is another regularization parameter. So the 
prior distribution in formula (4) is the product of (11), 
(12) and (13): 

(14) p pr (x(j)) = pr7V(x(j)) pos(xO)) pr3D(x(j)). 

The measurements are taken into account in the form 
of the likelihood distribution 

(15) 

where S is the covariance matrix of the Gaussian noise 
vector e. c ia a normalization constant. 

In reality, X-radiation measurement noise is Poisson 
distributed. Formula 15 uses a Gaussian approximation 
to the Poisson distribution. 

All the pans are presented of the right hand side of (4). 
Next it is looked for the MAP estimate of x(j), that is, 
the image x(j) that gives the largest value for the 
posteriori distribution (4). This is equivalent to finding 
the image x(j) that minimizes the expression 

F(xW) = hm^ - Ax^T-z-i^U) - Ax&) 4- 
- log(prTV(a;k'>)) - log(pr3D 



(16) 

with the additional requirement that every pixel of x(j) 
is positive. 



Minimization of (16) is difficult since F is not a 
differentiable function due to the absolute values 
present in (11) and (13) and due to the sharp cutting in 
(12). The absolute values in (11) and (13) are replaced 
by 

1*1 « M*) = |log(cosh(/?i)), 
(17) 

where p>0, enabling the use of efficient gradient-based 
minimization routines. The positivity constraint is taken 
care of by an interior search minimization method. The 
results are shown in figure 9 in comparison with 
traditional tomosynthesis technique based on 
backprojection. 

In this description explained mathematic modelling and 
other image information processing are performed by 
software means 15 to accomplish three-dimensional 
information of the object 4. 

In the second preferred embodiment the invention is 
utilized in dental panoramic x-ray imaging- The object 4 
is typically teeth of the patient and the medical x-ray 
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device 5 arrrangement is dental panoramic x-ray device 
arrangement 

In the third preferred embodiment the invention is 
utilized in mammography, for example in FFDM (Full 
Field Digital Mammography). There the object 4 is 
breast of human being and the medical x-ray device 5 
arrangement is mammography device arrangement. 

In the fourth preferred embodiment the invention is 
utilized in surgical x-ray imaging. There the object 4 is 
the body part under surgery and the medical x-ray 
device 5 arrangement is surgical x-ray device 
arrangement. 

In the first, second, third and fourth preferred 
embodiments of the invention the basic method steps 
are same as mentioned with the flow chart in figure 6. 
The utilizing of the invention in the second, third and 

fourth preferred embodiment is similar to what is 
described with the first preferred embodiment of the 
invention and elsewhere in this application except 
different medical x-ray imaging applications and their 
differences because of different medical x-ray devices 
and different objects to be x-ray imaged. 

What makes this invention especially inventive is that 
the imaging geometry and/or the motion of the object is 
solved with the image reconstruction simultaneously or 
having different timings. This is based on that projection 
data itself contains sufficiently information about the 
imaging geometry and also sufficiently information 
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about the 3D reconstruction. Neither reference feature 
points nor prior knowledge of the motion are needed. 

This can be done as described above with the preferred 
embodiments and with the means 15 for modelling the 
object 4 mathematically utilizing the projection data to 
solve the imaging geometry and/or the motion of the 
object and with the means 15 for utilizing said 
projection data and said mathematical modelling of the 
object in Bayesian inversion based on Bayes' formula 

P <. x ,e\m)= p " mp " ix)p{m \ x ' e) 

p(rri) 

to produce three-dimensional information of the object, 
the prior distribution p pr (8) representing the prior 
knowledge of the imaging geometry anci/or the motion 
of the object 4, the prior distribution p pr (x) representing 
mathematical modelling of the object, x representing 
the object image vector, which comprises values of the 
X-ray attenuation coefficient inside the object, e 
representing the parameter vector of the imaging 
geometry and/or the motion of the object, m 
representing projection data, the likelihood distribution 
p(m|x,G ) representing the X-radration attenuation 
model between the object image vector x and 
projection data m, geometry parameter vector 9, p(m) 
being a normalization constant and the posteriori 
distribution p(x,Q|m) representing the three- 
dimensional information of the object 4 and the 
imaging geometry including the motion of the object. 
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Although the Invention is described above with 
reference to the examples illustrated in the attached 
figures, it is to be understood that the invention is not 
limited thereto but can be varied in many ways within 
the inventive idea disclosed in the attached claims. 
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